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The locally self-consistent real space multiple scattering technique has 
been applied to calculate the electronic structure and chemical binding for 
the c(2x2)O/Cu(001) system, as a function of do-Cui ~ the height of oxy- 
gen above the fourfold hollow sites. It is found that the chemical binding 
between oxygen and copper has a mixed ionic-covalent character for all plau- 
sible values of do-Cui- Furthermore, the electron charge transfer from Cu to 
O depends strongly on do-Cui and is traced to the variation of the long-range 
electrostatic part of the potential. A competition between the hybridization of 
Cul-dxz with 0-px,Py and Cul-(i^2_j^2 with 0-pz states controls modification 
of the electronic structure when oxygen atoms approach the Cu(OOl) surface. 
The anisotropy of the oxygen valence electron charge density is found to be 
strongly and non-monotonically dependent on do-Cui- 

I. INTRODUCTION 

Oxygen adsorbtion on Cu(OOl) has been the subject of much discussion and debate in 
the past three decades mainly because of the illusive nature of the adsorption geometry. 
The history of the field is nicely summarized in a recent paper [jl]]. Here we mention some 
of the essential points. Unlike the case of Ni(OOl) in which a true c(2 x 2) phase is observed 
as the oxygen coverage approaches 0.5ML [0 in which the oxygen atoms occupy fourfold 
hollow sites at a distance about 0.8A above the top Ni layer, the stable surface struc- 
ture for similar oxygen coverage on Cu(OOl) has the periodicity (2^2 x V2)i?45° § and 
involves missing rows of Cu atoms according to LEED [^] and surface X— ray diffraction 
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0] measurements. This conclusion was agreed upon after much debate and analysis of data 
from several complimentary experimental techniques, some of which indicated the presence 
of a c(2 X 2) structure on an unreconstructed Cu(OOl) [[7|-[T0|]. Other experiments have in- 
dicated the presence of both of these phases |TT] , p^ and still others have given evidences of 
disordered phases in the system |]T3[. A more recent scanning tunneling microscopic (STM) 



study has also affirmed the presence of c(2 x 2) overlayer of adsorbed oxygen on Cu(OOl) at 
low coverage, but only in nanometer size domains and not large well-ordered domains . 



Currently, based on the above findings, the following picture can be built. If the adsorbate 
coverage is low, oxygen atoms adsorb on Cu(OOl) forming c(2 x 2) islands. Increase in the 
coverage leads to the ordering of the c(2 x 2) domains, forming c(4 x 6) pattern. When the 
coverage approaches a critical value of 34% of monolayer, the (2\/2 x \/2)i?45° structure 
made by a missing row type reconstruction is formed. 

The competition between the observed O/Cu(001) phases and the subsequent reconstruc- 
tion of the surface with the {2\^x ^/2)-R45° overlayer shows signs of a delicate balance in the 
surface electronic structure of the system which tilts in favor of one phase or the other de- 
pending on the coverage. A detailed description of the electronic states and chemical bonds 
formed when oxygen atoms adsorb on the Cu surface should thus provide much needed 
insight into the factors responsible for the observed variations in the surface structure. 

A possible theoretical explanation follows from effective medium theory fl^ which pre- 
dicts a reconstructed stable state for the substrate. These researchers find that the O valence 
levels interact very strongly with the Cu rf-bands such that the states can be shifted through 
the Fermi level. The difference in the total energy of the reconstructed and the unrecon- 
structed Cu surface arises from dependence of the energetic location of the antibonding 
levels, derived from the O 2p - Cu 3d hybridization, on coordination. A missing row can 
force this level up, thereby pushing up the hybridized antibonding level which can straggle 
the Fermi level, leading to a reduced occupation and thus stronger net bonding. In the case 
of a Ni substrate the 3d band straggles the Fermi level, and so its position is largely pinned 
by the need for charge neutrality. 



The rationale for a reconstructed Cu surface with O as adsorbate has also been related to 
the instability of the unreconstructed surface due to the large charge transfer to the adsorbate 



T6|,p^. Calculations performed for oxygen on small clusters of Cu (2-25 atoms ||T8| and 5 



atoms |T^) find a strong electronic charge transfer from the metal atoms to oxygen in the 
amount of 1.2e or 1.5e [|T^]. These two factors (the degree of the pO - dCu hybridization 
and amount of the effective charge on O) have thus far been proposed as the driving forces 
for the phase instability of the c(2 x 2) structure. The ab initio cluster calculations has 
also found the 0-Cu bond to be mostly ionic with some participation of the d-electrons of 



Cu ||T3- It should be borne in mind, however that models like the effective medium theory 
are at best semi-empirical. Further none of the theoretical methods as far employed take 
into account the contribution of the Madelung potential for these two overlayer structures, 
which may significantly change the subband energetic positions and charge distribution. If 
indeed an appreciable charge transfer and details of pO - dCu hybridization are responsible 
for the stability of a superstructure on Cu(OOl), a systematic calculation of these quantities 
is warranted. Such a calculation must also provide a comprehensive description of the 
relationship between the electronic and geometric structures of the system, including the 
contribution of the Madelung potentials. This paper is an effort in this direction. 

We have carried out detailed calculations of the electronic structure of the O/Cu(001) 
system using a reliable method and by varying the separation between the oxygen overlayer 
and the Cu(OOl) surface. As we shall see, there continues to be a debate also on the height 
of the oxygen atoms above the Cu surface. Our calculational strategy allows us to compare 
the local charge redistribution, modifications of hybridization of various electronic states 
and variations of the local and long-range contributions to the potentials, as a function of 
the height of the oxygen overlayer above the surface. It is through such systematic studies 
and, in conjunction with experimental observations, that we expect to gain insights into the 
nature of the chemical binding between the O and Cu atoms and establish rationale for the 
formation of a particular overlayer. Since, both c(2 x 2) and (2^/2 x ■\/2)-R45° phases have 
been shown to exist [|I|,|14|, our focus in this paper is on the former. At a later stage we will 



undertake a study of the missing row reconstructed surface with the large surface unit cell 
to represent the (2v^ x \^)RA5°, which would require a much larger computational effort 
than the present one. It should be mentioned that Wiell et al. [|T9[ have used a realistic 
full-potential LMTO method to examine the hybridization of orbitals of the N overlayer on 
Cu(OOl). These authors have also touched upon the O/Cu(001) system, but no details of 
the O - Cu electronic state hybridization or the charge transfer have been provided. 

The structural studies of the c(2 x 2) O/Cu(001) system performed with different tech- 
niques provide quite different results as the interlayer 0-Cu spacing (do-cui) ranges from 
to 1.5A. However, a majority of results [pl,P,p!2| gives do-cui value in the range 0.4- o.sA. 



The SEXAFS measurements ||T^ suggest three different configurations. In one of them the 
O atoms are located at the fourfold hollow (FFH) sites with do-cui = 

0.5A. Two other 

configurations represent the distorted FFH geometries with do-cui = 0.8A. The authors 
of recent photo-electron diffraction studies also propose a two-side model involving the O 
atoms at FFH positions with do-cui = 0.41 and 0.7A [|T|. They associate these two values 
with the edge and center positions in small c(2 x 2) domains, respectively. Our calculations 
are performed for the exact FFH geometry for a set of do-cui including the values obtained 
from the different experiments. We use a computational method based on the theory of 
multiple scattering in real space. This allows us to explicitly obtain the local electronic 
structure of atoms with different locations. Our results show a general picture of the elec- 
tronic state modification when oxygen atoms approach the Cu(OOl) surface and reveal the 
factors controlling charge transfer and chemical bond formation. 



II. COMPUTATIONAL METHOD 



The calculations are based on the density functional theory within the local density ap- 
proximation (LDA) 1^ and multiple scattering theory. Our approach embodies the local 
self-consistent multiple scattering method |^ , in which a compound is divided into overlap- 
ping clusters called local interaction zones (LIZs) centered around atoms in different local 
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environments. For each LIZ, a system of equations for the T— scattering matrix in the 
lattice site - angular moment representation for the muffin-tin (MT) potential |^2[ is solved 
self-consistently. The solutions are used to determine the cluster Green's functions and 
subsequently the local and total electronic densities of states and valence charge densities. 
An application of this method to materials related to the work here (copper oxides) has 
demonstrated its high efficiency and reliability pSj^ . 

In our model system, apart from the oxygen overlayer, the top four Cu layers (Cul - 
Cu4) were considered to be different from the bulk Cu. The fifth and lower Cu layers were 
supposed to have bulk properties. LIZ's were built around nonequivalent atoms belonging to 
the five different layers and contained 71 to 79 atoms depending on the local configuration. 
The sizes of the LIZ's were taken such that the calculated characteristics of the bulk system 
closely matched those obtained from other reliable methods. 

At each iteration of the self-consistent process a new potential is built by solving the 
Poisson's equation with the charge density obtained at the previous iteration. To take into 
account surface effects we follow the approach developed in Ref. where the MT potential 
at the ath atom belonging to the ith layer is given by 
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(1) 



The nuclear and the MT electronic charge density contributions to the potential (the first 
term and the terms in square brackets in (1), respectively) are of the same form as in the 
bulk calculations p6|j27|] . The electrostatic (V^'-'') and exchange-correlation (V^^^(p°)) parts 



of the MT constant are taken from the bulk. The last two terms in the above expression 
represent the distant part of the MT potential: 



Vd = -Anp^ + V^J'^'' (2) 

This part contains explicit information about the surface. The interstitial charge density 
p° is supposed to be layer-dependent reflecting the asymmetry of the system. In such an 



approach the space in the vicinity of the surface is divided by layers belonging to different 
atomic planes parallel to the surface. In the present calculations we keep the thickness of 
the oxygen layer independent of the interlayer O - Cul spacing such that the radius of 
the corresponding Wigner-Seitz sphere is equal to the O-II ionic radius Rq = 1.24A. The 
thickness of the Cu2 and next Cu layers corresponds to the bulk lattice parameter. The Cul 
layer takes the rest of the space between the O and Cu2 layers. The p° values are obtained by 
averaging the interstitial charges. The latter are differences between Wigner-Seitz (Q^^) and 
MT {q^a'^) charges, which are calculated by integration of the electronic charge density over 
Wigner-Seitz and MT spheres, respectively. The values are multiplied by a coefficient 
to meet charge neutrality condition for the region disturbed by the surface. This coefficient 
approaches unity, with charge convergence. The Madelung potential includes monopole 
C,.^, dipole Mll^^ 



M^^ja, dipole and interstitial Vij[pj] terms |25 



V^' = E (?J73^S.73 + d,,Mll,, + V.J [Pj]) , (3) 



where dji3 is a dipole moment of the MT charge density and 
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The M^^jp and Ml^j^j terms in (3) are calculated by an Ewald-like plane-wise summation, 
and explicit solution of one-dimensional (normal to the surface) Poisson's equation, using 



the technique described in Ref. p5| 



The exchange and correlation parts of the potential are determined within LDA using 



the technique described in Ref. |28 



To test the applicability of the method to surface problems, we have calculated the 
electronic structure of the clean Cu(OOl) surface and compared our results with ones obtained 



by means of TB-LMTO method p9|, screened KKR Green's function method BO] and self 



consistent localized KKR scheme developed for the surface calculations ||2^ . As we shall see 
in the next section, our projected densities of electronic states calculated for the c(2 x 2) 
O/Cu(001) system are in good agreement with those obtained by full potential LMTO 
method [|l9l 



III. RESULTS AND DISCUSSION 



We have calculated the electronic structure of the O/Cu(001) surface with the c(2x2) 
symmetry for selected values of the 0-Cul interlayer spacing do-Cui- As already dis- 
cussed, there is some discrepancy in the values of do-Cui obtained from different experi- 
ments. However a majority of experiments have found do-cui to lie between OAA and O.SA 
(0.11a — 0.22a, where a is the lattice parameter of bulk fee copper). With this in mind we 
have varied do-Cui from 0.15a to 0.3a in our calculations and examined its effect on the 
surface electronic structure. The local densities -A^"(-E') of electronic states, their projection 
A^/^(i?) on the cubic harmonics, and the valence electron density p{r) have been calculated 
for the oxygen overlayer and the four top Cu layers, for selected values of do-Cui- 

The densities of the pO- and (iCu-electronic states of the system are plotted in Fig. 1 
for the longest and shortest interlayer spacings considered. As seen in Fig. 1, the N^'^(E) 
calculated for the Cu3 and Cu4 atoms (third and fourth layer) are hardly affected by the 
choice of do-Cui and differ slightly from the bulk values whereas the (iCu2-states are notice- 
ably disturbed. The dCul and pO state densities are found to have a rich structure and for 
both values of do-Cui the spectra show a clear signature of the pO - dCul hybridization. 
Note the position of the peaks marked "a" and "6" for do-cui — 0.3a and "a"' and "6"' 
for do-cui = 0.15a in Fig. 1. When the oxygen overlayer approaches the Cu(OOl) surface 
{do-Cui = 0.15a), the splitting (the distance between a' - b' peaks in the figure) increases 
reflecting an enhancement of the hybridization, which even involves the second copper layer. 

To examine the charge transfer in the vicinity of the surface we have calculated the 
Wigner-Seitz charges that allowed us to obtain the amount of charge deviation (AQj) from 
electric neutrality per 2D— unit cells belonging to the different layers. The results are shown 
in Fig. 2. One can see a strong charge transfer from the Cul layer to the O overlayer 
coupled with some increase in the electron charge in the Cu2 layer. The effect is significantly 
enhanced when the oxygen atoms approach the Cul layer. From this result one would 
conclude that the O - Cul chemical binding has an essentially ionic character and that the 
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ionicity increases as oxygen atoms approach the Cu(OOl) surface. In Table 1, we present 
the amount of charge transfer for the five different heights of the oxygen overlayer that we 
considered. At the shortest 0-Cul distance, the two Cu atoms in the unit cell of the first 
layer have a deficit of about one electron each, with 1.5 electron going to the O atom and 
0.25 electron for each of the 2 atoms of the second layer. The amount of charge provided to 
the O atom drops to 1 electron when it's height above Cul is doubled to 0.3a. The charge 
deviation AQo of the O layer and the distant part Vd of the potential in Eq. (2) at the O 
site are plotted versus do-Cui in the lower panel of Fig. 2. A clear correlation between AQo 
and Vd seen in the figure suggests that the long-range electrostatic interaction is a driving 
force for the charge transfer in the system. 

In the system under consideration, each oxygen atom and its four nearest Cul neighbors 
form a pyramid. Such an atomic configuration is expected to cause a strong anisotropy in 
the electronic structure. Therefore an analysis of the N^{E) distributions for O and Cul is 
very useful for understanding the nature of the chemical bonds formed when oxygen adsorbs 
on the Cu surface. A simple symmetry consideration shows that only the 0-PxiPy 

-Cul-d,, 

and 0-pz - Cul-dx2-.y2 hybridization can be significant in the O - Cul subsystem. Therefore 
our focus is mainly on these electronic states. In Fig. 3, the densities of the 0-pz and Cul- 
dx'2-y'2 states are plotted for the values of do-cui as labeled. Both projected densities are 
considerably modified with the decrease in do-cui- The band broadens and distinct new 
peaks appear as do-Cui takes the values 0.2a and 0.15a. An energetic alignment of the 
0-pz and Cul-dx2_y2 peaks in the a, b and c regions is a signature of hybridization of these 
states. An increase in the peak amplitudes and an extra splitting with the decrease in do-Cui 
reflect a strong enhancement of the hybridization as the oxygen atoms approach the surface. 
In contrast, the splitting of the 0-px,Py and Cul-d^z states is almost independent of the 
do-Cui value (see Fig. 4) . This means that the hybridization of these states is not changed 
noticeably, when do-Cui is varied within the range that we considered. Such a difference 
in the response of the electronic states to variation of the 0-Cul interlayer spacing can be 
explained on the basis of a simple geometric analysis. When the oxygen atoms approach the 
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Cul layer, the O - Cul bond length decreases and the symmetry of the 0-pz and CvLl-d^2_y2 
cubic harmonics is such that their overlap increases enhancing the hybridization. On the 
other hand, the symmetry of the 0-px,Py and Cul-d^z states is such that the overlap of 
their cubic harmonics is diminished with decrease in do-cui that compensates the effect of 
the O - Cul bond length reduction. 

Thus, according to these results, a noticeable covalence pO - dCul binding takes place 
in the system along with the strong ionic binding. If the O - Cul interlayer spacing is long 
enough {do-cui = 0.3a), the covalence component of chemical binding is mostly determined 
by the 0-px,Py - Cul-dxz hybridization. When do-cui is smaller, the 0-pz - Cul-dx2-y2 
contribution to the covalence bond increases. This competition mostly controls modification 
of the electronic states when the oxygen atoms approach the surface. 

The densities of 0-px,Py and 0-pz electronic states have also been calculated for 
O/Cu(001) by means of the full-potential LMTO method [0, which is considered to be 



one of the best approaches for electronic structure calculations. These authors considered 
an oxygen overlayer with a c(2 x 2) structure with the oxygen atom occupying an FFH and 
a quasi-FFH sites. In Fig. 5 we have compared our results obtained for do-cui = 0.15a 
(0.514A) with the results of Ref. |I§ at the FFH site at a height of 0.5A, for 0-pz and 



0-Px,y Both the splitting and peak amplitudes are in very good agreement attesting to the 
reliability of the method used here. We note from Fig. 5 that a shift the height between 
0.15a and 0.2a produces a shift of about leV of the spectrum, again in good agreement with 
the results in Ref. |^ . 



We have calculated the radial dependence of the valence electron charge density around 
the O and Cul - Cu4 atoms along some symmetrical directions. The plot of the oxygen 
charge density in Fig. 6 displays it to be essentially anisotropic. We have also found 
that the anisotropy increases when do-cui is reduced from 0.3a to 0.2a whereas further 
do-Cui reduction leads to a decrease in the anisotropy. This effect is illustrated in Fig. 7, 
in which the radial dependence of the difference between the densities along the (001) and 
(100) directions is plotted for systems with different values of do~cui- Such a non-monotonic 



behavior can be explained by the hybridization competition mentioned above. The change in 
the anisotropy of the oxygen charge density occurs because different amounts of charge flow 
to the 0-px and 0-pz states, as do-cui decreases. The stronger covalence binding, the larger 
an increase in the electron-electron repulsion caused by an extra electronic charge in the 
bond region. Therefore the charge tends to come to the states less involved in hybridization, 
which are the 0-pz states when do-cui is large. This leads to an increase in the anisotropy. 
The tendency is changed at short interlayer distances, when the hybridization of the 0-pz 
and Cul-dx'i-yi is enhanced significantly that can make preferable an occupation of the 
0-Px,Py states. 

We now turn to the implication of the results presented here to a simple model proposed 
for the adsorption of oxygen on Cu(OOl). To describe this phenomenon, the authors of Ref. 



T4 use a model isotropic potential based on the local part of the Madelung potential. At the 
same time they assume that not only the O and Cul layers are involved in charge transfer, 
but also deeper ones. Our results confirm the latter assumption (see Fig. 2) and imply that 
an accurate model potential should be based on the long range Coulomb interaction rather 
than on the local one. Furthermore, using the experimental data obtained for the c(2 x 2) 
and (2-\/2 x \/2)RA5° phases the authors of Ref. |]T^ build the isotropic potential, which 
appears to have step-like dependence on the O - Cul bond length. To obtain this, these 
authors attribute a step-like behavior to the effective atomic charges. They assume that 
such a behavior can be induced by an increase in covalence, as the O - Cul bond length is 
reduced. Our results do not support these assumptions. We have found that the charges and 
covalence of the bonds are smooth functions of the interatomic distance. The c(2 x 2) and 
(2\/2xy2)i?45° phases are different in the symmetry of the oxygen local surrounding rather 
than in the interatomic distances. This suggests that the phase transition can be properly 
described only by means of an anisotropic potential. Moreover, a noticeable covalence of the 
O - Cul binding (even for relatively long bond lengths) and its anisotropy obtained in our 
study allows us to propose that an accurate modelling of the c(2 x 2) phase itself requires 
an anisotropic potential. 
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IV. CONCLUSIONS 



The electronic structure has been calculated for the O/Cu(001) system in the c(2 x 2) with 
five different O - Cul interlayer spacings. We find that the oxygen and copper atoms form 
a mixed ionic-covalence chemical binding in the O/Cu(001) system for the whole range of 
the considered O - Cul interlayer spacings. The oxygen atom gains an extra charge varying 
between 1 and 1.5 electron when its height is varied between 0.3a and 0.15a. This charge 
comes solely from the atoms at the top layer of the substrate with a small enhancement of 
the second layer electronic charge. We also find that this charge transfer is controlled by 
variation of the long-range electrostatic interaction. The detailed analysis of the properties 
of the local electronic structure reveals that the electronic structure of the O/Cu(001) system 
is governed by a competition between the hybridization of Cul-d^z - O-p^-iPy a-nd Cul-<i-r2_j,2 
- 0-pz states, which depends on O - Cul spacing. The anisotropy of the oxygen valence 
electron charge density is found to be strongly and non-monotonically dependent on the 
interlayer spacing. 
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FIG. 1. The density of p-electronic states of the oxygen over layer and d-electronic states of the 
four top Cu layers for O/Cu(001) for 0-Cul interlayer spacings of 0.3a (dashed line) and 0.15a 
(solid line). 

FIG. 2. Upper panel: the layer charges per 2D-unit cell. The layer ^ corresponds to the 
oxygen overlayer. Lower panel: charge deviation for the oxygen overlayer and the distant part of 
the MT potential plotted versus the 0-Cul interlayer spacing. 

FIG. 3. The densities of the 0-pz (dashed line) and Cul-d^2_y2 (solid line) electronic states for 
c(2 X 2) O/Cu(001) with different 0-Cul interlayer spacings. 

FIG. 4. The densities of the 0-px,y (dashed line) and Cul-da;^ (solid line) electronic states for 
c(2 X 2) O/Cu(001) with different 0-Cul interlayer spacings. 

FIG. 5. The densities of the 0-px,y and 0-pz electronic states calculated in the present work for 
two O heights (solid line 0.541 A) and (long-dash dotted line 0.722 A) compared with those in Ref. 
18 (dashed line) for c(2 x 2) O/Cu(001) at 0.5 A. A convalution with a lorentzian with 0.7 eV full 
width at its half maximum is used. 

FIG. 6. The oxygen valence electron charge densities plotted along the (001) (dashed line) and 
(100) (solid line) directions for O/Cu(001) with two different 0-Cul interlayer spacings. 

FIG. 7. The radial dependence of the difference between the oxygen valence electron charge 
densities calculated along the (001) and (100) directions for O/Cu(001) with three different 0-Cul 
interlayer spacings. 

TABLE I. Valence electronic charge variation A.Q{e) for different 0-Cul interlayer spacing. The 
percentage change is given in parentheses 

Layer 0.15a 0.17a 0.20a 0.25a 0.30a 

O 1.493 (37.3%) 1.429 (35.7%) 1.329 (33.2%) 1.144 (28.6%) 0.983 (24.6%) 
Cul -2.108 (-9.6%) -1.986 (-9.0%) -1.800 (-8.2%) -1.484 (-6.8%) -1.211 (-5.5%) 
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Cu2 0.575 (2.6%) 0.519 (2.4%) 0.430 (2.0%) 0.313 (1.4%) 0.207 (0.9%) 
Cu3 0.040 (0.2%) 0.039 (0.2%) 0.041 (0.2%) 0.027 (0.1%) 0.021 (0.1%) 
Cu4 0.000 (0.0%) 0.000 (0.0%) 0.000 (0.0%) 0.000 (0.0%) 0.000 (0.0%) 
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